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The dark energy in the Universe is described in the context of modified Friedmann equations as 
a fluid parameterized by the density of dark matter and undergoing an adiabatic expansion. This 
formulation is applied to the Cardassian model. Choosing then parameters consistent with the su- 
pernova observations, it gives a background expansion in which the cosmic temperature fluctuations 
are calculated. The resulting spectrum is quite similar to what is obtained in the standard concor- 
dance model. If the Cardassian fluid is interpreted as a new kind of interacting dark matter, its 
overdensities are driven into oscillations when the interaction energy is rising in importance. This 
does not occur in a description of the Cardassian fluctuations motivated by theories of modified 
gravity. There the energy of the underlying matter is also conserved, which requires appearance of 
effective shear stress in the late Universe. In both approaches that allow fluctuations the thermal 
power spectrum at large scales is much too strongly enhanced by the late integrated Sachs- Wolfe 
effect. With the interacting dark matter assumption, we conclude that the Cardassian model is 
ruled out by observations, expect in a small neighbourhood of the ACDM limit. 



Introduction 



parameter a(t) in a spatially flat Universe, 



Observations of distant supernovae indicate that 
the expansion of the Universe now undergoes an 
acceleration^, 0, Q- This is also consistent with the 
most recent and accurate measurements of the fluctua- 
tions in the cosmic microwave background by the WMAP 
satellite0|. These also imply that the Universe is at least 
very nearly spatially flat. 

A possible explanation of these observations requires 
that the Universe contains a new and unknown com- 
ponent, called dark energy, in addition to radiation, 
baryons and dark matter. This could be Einstein's 
cosmological constant or a more dynamical component 
based on the universal presence of a scalar fleld called 
quintessence^, 6J. Alternatively, one could contemplate 
the possibility that standard Einstein gravity is not valid 
at very large scales so as to allow for modiflcations of the 
Friedmann equations 0, 

We will here consider one particular such proposal by 
Freese and Lewis where the total energy density is 
made up of just the usual radiation, baryons and cold 
dark matter densities pr, pt and pc but they enter in a 
non-linear way in the Friedmann equation for the scale 
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where H — a/a is the Hubble parameter. There 
are several such proposals based on ideas from brane 
phvsicsflfij or some unknown interactions between mat- 
ter particles [rH- But in order to make contact with the 
CMB physics, we will deflne the model so that the Fried- 
mann equation takes the form 
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where M = (SttG) 2 is the reduced Planck mass and 

PK = Pc[l + Bp-'^'']'/^ (3) 

is the modified polytropic Cardassian energy density. 
Here _B is a positive constant, q is the polytropic index 
and it is modified as long as the parameter u ^ 1. All 
the unknown physics is now lumped into the dark sector. 
At early times when the matter density pc is large, the 
last term in Q will be negligible and px ~ Pc- How- 
ever, at late times when pc becomes small, the last term 
dominates and will act as a dark energy component due 
to unknown properties of dark matter. The energy frac- 
tions fli = Pi/ Per with Per ~ iA'PH^, are seen from (jSJ 
to satisfy 51^ + 5^6 + = 1 since there is no curvature 
in space. In our calculations we use the WMAP valuesQ 
for these parameters. Then Oc = 0.226, fib = 0.044, and 
fir we get from assuming three massless neutrino species 



2 



and the background temperature T = 2.726 K. For h we 
use the value 0.72. Today we will therefore have flKo — 1 
when we neglect the small contributions from radiation 
and baryons. Given q and v, the remaining parameter B 
follows then from 



(4) 



Ordinary dark matter has zero pressure and satisfies the 
conservation equation pc + 3Hpc = (in the homoge- 
neous and isotropic background Universe) . Thus we have 
Pc — E^nd correspondingly for the other standard 

ingredients, with a — 1 today. The Friedmann equation 
(O can then be integrated to give the full background 
evolution of the scale parameter a{t) with Cardassian 
energy present. 

Following Gondolo and Freese|lJ, we will assume that 
the Cardassian energy density px is that of a fluid which 
undergoes an adiabatic expansion when the Universe 
evolves. It will thus have a pressure pk which can be 
obtained from its conservation equation pK + 3H{pK + 
Pk) = 0. Eliminating the Hubble parameter H between 
this and the conservation equation for the dark matter, 
the Cardassian pressure will follow from 



dpK 

PK = Pc—, pK- 

dpc 

For the choice (|2Jl it gives 

PK = -i^BpZ'^''+^ \l + Sp-'"- 
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This equation of state can be written on the standard 
form Pk — wkPk with 
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At late times when pc — > we see that wk —v. We 
therefore expect late-time acceleration when the parame- 
ter V > 1/3. The equation-of-state parameter wk is zero 
at early times, and approaches the value wk = —v in the 
late universe, more rapidly for larger q. Late evolution 
of Wk for different values of the parameters q and v is 
shown in FigQ] 

The late expansion history of these models has been 
compared to the observational data of su p ernovae [l^ 
and of other astrophysical objects The locations 
of acoustic peaks in the CMB spectrum has also been 
considered before0,0|. We wiU calculate the full CMB 
spectrum in these models, first assuming that the Cardas- 
sian component accounts only for the background expan- 
sion of the Universe. This would, however, correspond 
to a time-variable but nonfluctuating vacuum energy, to 
which we have no physical motivation|34|. This is our 
case I. 

A more realistic description, which allows fluctuations 
in the Cardassian fluid, is our case II. Then the Car- 
dassian fluid can be thought of as an interacting dark 
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FIG. 1: Equation-of-state parameter wk as a function of the 
scale factor. The upper panel has v = \ and g-values q = 
1,2,3,4,5 and g = 6 (the steepest curve is for g = 6). The 
lower panel has g = 1 and v — 0.6,0.8,1.0,1.2, 1.4 and v — 1.6 
(the most negative wk today is for v — 1.6). The dashed 
curves correspond to ACDM. 



matter. The energy densityj35| pK — Pc would then be 
attributed to an unknown particle species or field me- 
diating the dark matter interactions. Thus Pc does not 
satisfy an independent energy conservation law. It only 
satisfies a mass conservation law, whereas we have energy 
conservation for the total pk- 

In our case III we impose energy conservation sepa- 
rately for the cold dark matter. This can be motivated 
by theories of modified gravity. In these theories, we can 
generally write the Einstein equations as equations as 



n — ^ (rrM 
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where M is for matter and C for the corrections to the 
standard gravity. The latter can be parameterized as 
a function of the matter energy density. We illustrate 
this with the example of corrections in the form of DGP 
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This can rewritten as 
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The RHS we then interpret as an energy density of a 
fluid. Granted that the matter conservation also now 
apphes as usually, we can find the pressure of this fluid 
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which gives the equation of state for this effective mat- 
ter. It would also follow by using Eq.jSJ with Ea. (|10(l . 
However, in this paper we consider only Cardassian mod- 
ifications to the Friedmann equation. For a recent inves- 
tigation of the cosmological constraints for modifications 
generalized from Eq.®, see Elgar0y and Multamaki[T8l|. 
Note also that wc have included only the contribution 
from the density of dark matter in Eq.Q, because the 
corrections to the Einstein gravity do not become impor- 
tant until the matter dominated era. This justifies our 
use of Eq.(|21l also in the case IIl|3^. 



Case I: No Cardassian fluctuations. 

In order to calculate the temperature fluctuations in 
this Universe, we now consider these cases. Firstly we 
assume that the Cardassian fluid just provides a modified 
background for the evolution without any internal fluc- 
tuations. These take place in the ordinary components of 
radiation, baryons and dark matter and can be calculated 
by standard methods for any values of q,!^ > 0. In FigEl 
we show the Bardeen potentials for the choice v = 0.8 
and q = 1.5 and for different values of the wave number 
k which sets the scale of the fluctuations. [s^ The evolu- 
tion of gravitational potentials is seen to be very similar 
to what is found in ACDM models on all scales. 

We have calculated the CMB spectra using the fluid 
approximation for photons and the well known analyti- 
cal result for the diffusion damping scale (e.g. [IqS). We 
have tested the fluid approximation with various models 
and found it to be in agreement with more exact calcu- 
lations (using e.g. CMBFAST"20'|) within 5 per cent for 
£ smaller than a few hundred. This is sufficiently accu- 
rate to uncover the interesting large scale features of the 
models we are considering, as is the main purpose of this 
investigation. For simplicity, we do not include reion- 
ization. Our perturbations are normalized such that the 
primordial curvature perturbation TZ = 1. 




FIG. 2: Evolution of * (solid line) and $ (dashed line) in 
a Cardassian background (case I). From bottom to top, k — 
0.00024, 0.0012, 0.0024, 0.0060, 0.012, 0.036, 0.048 Mpc-\ 




FIG. 3; CMB spectra in a Cardassian background (case I). 
From bottom to up, (q, i/) = (1, 1), (1.5, 0.8), (2, 0.7), (10, 0.6) 
and (100, 0.6). The bottom curve with both q and u equal to 
1 corresponds to the ACDM model. 



In Fig|31 we show the full thermal spectrum 
obtained from this modified background. Here 
we use the same Cardassian parameters (q, v) = 
(1,1), (1.5, 0.8), (2, 0.7), (10, 0.6) and (100,0.6) which 
have been found by Savage, Sugiyama and Freesefl^ to 
be consistent with the SNIa observations and the age of 
the Universe. The peaks in the thermal spectra are lo- 
cated at the right places and the contribution from the 
integrated Sachs- Wolfe effect is seen to increase modestly 
with q. However, by lowering q below 1 one cannot re- 
duce the ISW effect without losing the agreement of the 
spectrum with WMAP observations on smaller scales. 
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Case II: CDM fluctuations driven by the Cardassian 

fluid. 



As our second approach, we include fluctuations in 
the Cardassian fluid. These wiU then drive the fluc- 
tuations in the underlying dark matter. In our dis- 
cussion of perturbations we adopt the notation of Ma 
and Bertschinger(22j . with the only exception of using 
uppercase letters for the gravitational potentials in the 
conformal-Newtonian gauge. In this gauge the line ele- 
ment is then written as 



a2(r)[-(l + 2*)dr2 + (1 - 2'P)S,jdx'dx^]. (13) 



We assume adiabatic perturbations, which means that 
the relations © and © between pK, PK, and pc will 
hold also for the perturbed quantities. 

Assuming no anisotropic stress in the Cardassian fluid, 
we find that its density perturbations in this gauge obey 
the equation of motion 

S'K = il+ WK){-eK + 3$') + 3n{wK - cI)Sk (14) 

while the corresponding velocity perturbation is governed 
by 



{3wK - i)neK 



1 + Wk 



1 + Wk 



(15) 



Here the prime denotes the derivative with respect to 
conformal time, H, — a' /a and the equation of state pa- 
rameter Wk is given by Q . Its conformal time derivative 
is then 



iUqv'^Bp-i'' 
" [1 + Bp-^y 



(16) 



and goes to zero at late times. The remaining parameter 
is the Cardassian speed of sound 38] c\ = {dpK/dpK)s = 
dpK/dpK- From the pressure (jSJ wc then find 



uBp-'i'^[{v-l)Bp, 



qu 



qv 



-K 



1 + (2 - v)BpT + (1 - v)B-^p, 



(17) 



At late times when ^ this is seen to approach c|. — > 
—V which is unacceptable. This can be avoided only by 
choosing the special value v — \. Then the expression 
for the sound speed simplifies to 



g-1 
1 + 



(18) 



In order for it to be positive, we must have 9 > 1 and 
also g < 2 so that c^^ < 1 at late times. 

We find it also useful to do the calculations in the 
synchronous gauge, in which the line element is given by 



ds^ = a^(r)[-dr^ + {6ij + hij)dx''dx'] 



(19) 



where we define the scalar modes in the Fourier space as 



/iy(x, r) 



kikjh(k, t) 



1 



Q(hkj - -(5ij )77(k, r) 



(20) 



when k = kk. In this gauge the fluctuation equations for 
a general fluid, which may have also anisotropic stress cr, 
are 

5' ^-{l + w){e + ^) + m{w-c^)5 (21) 



for the density contrast and 

nn' 

e' = (3w - i)ne 



\ + w 1 + w 
for the velocity perturbation. 



k^S-k^a (22) 




FIG. 4; Evolution of * (solid line) and <1> (dotted line) in the 
presence of a Cardassian fluid (case II). From bottom to top, 
k = 0.00024, 0.0012, 0.0024, 0.0060, 0.012, 0.036, 0.048 Mpc~\ 



We have numerically solved the perturbation equations 
in both gauges now including also the fluctuations in the 
Cardassian fluid for the representative value g = 1.5. In 
this case the fluctuations in the dark matter are not inde- 
pendent. The evolution of the Bardeen potentials ^ and 
<& at different scales is shown in Fig^J They are seen to 
have a very rapid increase (decrease in absolute value) 
at late times for all scales comparable and smaller than 
the present horizon. This results in an unusually strong 
late-time ISW effect. In the thermal spectrum in FigO 
it is seen to give a rather big deviation from the WMAP 
data for scales £ < 100. We find the same results in both 
gauges. The top curve is an example of a model where 
the sound speed of the Cardassian fluid exceeds the speed 
of light. Models where c|. takes negative values are not 
computable at all since fluctuations then blow up, except 
for tiny deviations from the ACDM values. 
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FIG. 5: CMB spectra in the case II. Here u = 1 and from 
bottom to up, g = 1, 1.05, 1.1, 1.2, 1.4, 1.8 and 2.6. 
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FIG. 6: Evolution of Sk in the case II. The model 
is the same as in Fig. 4. From bottom to top, k = 
0.0012,0.0024,0.0060,0.012,0.036,0.048 Mpc'^ Dashed 
lines represent the perturbations in Newtonian and the solid 
lines in the synchronous gauge. 



The dynamics is now different from the case I, where 
the growth of matter perturbations was slowed down only 
because of the accelerated expansion of the background 
Universe. Now the overdensities, which have grown in 
the Cardassian fluid during its matter-like behaviour, are 
swiftly decaying as the fluid is turning into an effective 
cosmological constant. Meanwhile, the decaying over- 
densities perform rapid oscillations due to presence of 
the Cardassian pressure. This is shown in Fig|B| Fluc- 
tuations in the cold dark matter and in the Cardassian 
fluid are related adiabatically 

S^^(^J^y\,^^^a^.S,, (23) 
\ d log pc J l + WK 



where the last equality holds at late times (and when v = 
1). From this we see that the dark matter perturbations 
are driven along with the Cardassian oscillations, but are 
not sharing the decay rate of the host fluid perturbations. 

The Cardassian fluid which is restricted here to have 
the two parameters v = \ and 1 < g < 2, is related to 
the generalized Chaplygin gas (GCG),23]. In the latter 
model the energy density is 

PocG^[B + ^^]^^. (24) 

After the change of variable, g = 1 + a, we see that this 
reduces to the Cardassian model in the special case that 
;/ = 1 (A is fixed by requiring standard early cosmology). 
In the Cardassian model we consider cold dark matter ex- 
isting separately and that its perturbations obey Ea. (|23|l . 
But the resulting CMB spectrum should be equivalent 
regardless of this intrinsic decomposition of the fluid. In- 
deed, our results agree with^j, where the CMB spec- 
trum was calculated for the generalized Chaplygin gas. 
Note that there the normalization of the spectra is dif- 
ferent from ours, since there the Sachs- Wolfe plateau is 
kept low, resulting in lower amplitude at small scales. 

For completeness, we consider also the matter power 
spectra in these models. We define here the total matter 
power spectrum, including both the contribution from 
the Cardassian fluid and from the baryons, as 

(25) 

27r^ 27r"^ V PK + Pb J 

In the FigE|we see that modifying only the background 
evolution has little effect on the matter power spectrum. 
However, as shown also by Amarzguioui, Elgar0y and 
Multamaki l25ll . consistently with Sandvik, Tegmark and 
Zaldarriaga[2q. when one adopts the fluid interpretation 
of the Cardassian expansion, even tiny changes in the 
parameters q^v result in observable departure from the 
ACDM cosmology in the matter power spectrum. This 
is seen in Fig|Sl When c| < 0, the spectrum blows up. 
When c| > 0, the spectrum of Cardassian fluid would 
be oscillating as expected from Fig|Hl but these are not 
seen in the total matter power spectrum since 5k is sup- 
pressed and thus the baryonic fluctuations dominate the 
structure at smaller scales. Also the amplitude of the to- 
tal matter power spectrum is sensitive to the parameters 
q and v. There has been some interesting suggestions on 
how to render the growth of structure in the GCG model 
more plausible [27i. i28i]. The one proposed inj22j imposes 
intrinsic entropy perturbations in the fluid to cancel the 
effect of the finite sound speed. This could be also done 
to the Cardassian fluid we are considering. In practice it 
would amount to simply setting c|- = in Eas. ll4|l and 
(|15|l . However, since entropy perturbations are gauge- 
dependent, it would perhaps be difficult to find physical 
justification for this. Another recent proposal 28] intro- 
duces an unique decomposition of the Chaplygin gas into 
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FIG. 7: The total matter power spectra for different param- 
eter values in case I. 



FIG. 8: The total matter power spectra for different param- 
eter values in case II. 



interacting vacuum energy and cold dark matter. In the 
case of Cardassian fluid, the corresponding decomposi- 
tion would require the redefinition of the cold dark mat- 
ter density as 

Pc = Pk + Pk 

= pS + Bp-'^'^]'.-'[l + (l-i.)Bp-n (26) 

Then we would find that the equation of state of the fluid 
Pk — Pc would be a constant, with w = —1. In|28l| it was 
assumed that such a component would be unfluctuating 
and that the redefined cold dark matter would have no 
pressure perturbations. However, these assumptions are 
invalid if one is restricted to standard general relativity 
and adiabatic perturbations, since then p'^^c — p'c^k, 
and now p'j^ ^ because of the interaction. 

Both of these proposals, the imposement of entropy 
perturbations and the decomposition of the fluid into 
dark matter and unperturbed vacuum energy, would 
prevent oscillations in the overdensities and result in a 
smaller ISW effect for the Cardassian model than in 
Figl^l but admittedly these proposals seem rather ad hoc 



unless one could specify a physical reason responsible for 
neglecting the undesiderable features in the evolution of 
inhomogenities in the fluid. Now we instead move on to 
consider the possibility that the form of the Cardassian 
Friedmann equation arises not from energetics of physical 
fluids but from modifications to the Einstein gravity. 



Case III: Cardassian fluctuations induced by the 
CDM. 

To avoid the dark matter oscillations driven by the 
Cardassian pressure, one would have to modify the sce- 
nario so that the dark matter would not see the Car- 
dassian pressure. This can be achieved, without relaxing 
the usual mass conservation or the adiabaticity of the 
cold dark matter, if it would satisfy energy conservation 
separately. This would remove any but the gravitational 
interaction between p^ and p^ — Pc, and so there should 
be no reason in ordinary physics why the parametric re- 
lation between the two should be maintained. However, 
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we have tried to motivate such a scenario by modified 
gravity as discussed in the introduction. 

Now both the matter and the Cardassian perturba- 
tions would obey the corresponding Eqs.j^U and (|22|l . 
In our case I they were satisfied only for the dark mat- 
ter, in our case II only for the Cardassian fluid. It is not 
immediately clear that they can be satisfied for both com- 
ponents simultaneously. This is easiest to see in a syn- 
chronous gauge, since in such particular choice of gauge 
the velocity perturbations in matter, and thus also in the 
Cardassian fluid can be set to zero. Thus we have only 
to solve Ea. H21(l for dark matter, 



h/_ 



(27) 



and invert Ea. iP^ to get the density perturbation of 
the Cardassian fluid. To check that there indeed are no 
Cardassian velocity perturbations in this gauge, one can 
rewrite Ea. (|21|l as 



h/_ 

T 



1 -I- Wk 



1 -I- Wk 



(28) 



Since by Eq. 123) S'j^ = SHiwK - cj^)SK + il + WK)S'^, the 
RHS vanishes identically. Eq. (|22l) tells us that there now 
is anisotropic stress in the fluid, which is proportional 
to the density perturbation evaluated in the synchronous 
gauge: 



OK 



1 -I- Wk 

[1 + (1 - iy)Bp-'y 



Sk. (29) 



Except for the last equality, this formulation applies gen- 
erally for modifled Friedmann equations with pxipc), 
since it does not depend on the actual form of Eq.©. 

Now the overdensities in pressureless dark matter give 
arise to perturbations in the Cardassian density, in which 
the shear stress acts to eliminate the effect of pressure 
gradients (and thus the sensitivity to c^). Thus no late- 
time oscillations occur in the overdensities of either com- 
ponent. In this case one has the freedom to choose any 
positive value for q and v, since Sk is now given by 
Ea. (|23|l and thus can take any value, even negative. 
The ISW effect is again strong, since now the Cardas- 
sian fluctuations at late times are suppressed according 
to Eq.^. 

The stress term becomes important at late times, and 
it causes deviation between the gravitational potentials 
in the Newtonian gauge. This is seen from the constraint 
equation 



^ = 



2 APk^ 



{p+p)a 



(30) 



3 a^vBp-i''+^ {v - l)Bp-i'' 



qv — \ 



j(Syn) 




FIG. 9: Evolution of (solid lines) and $ (dashed lines) in 
the presence of the Cardassian fluid (case III). From bottom 
to top, k = 0.00024,0.0012,0.0024,0.0060,0.012,0.036,0.048 
Mpc-\ 




FIG. 10: CMB spectra in the case III. Here u = 1 and from 
bottom to up, g = 1, 1.05, 1.1, 1.2, 1.4, 1.8 and 2.6. 



The a in the first line is the total gauge-invariant shear 
stress which at late times gets a contribution only from 
gk- In the second line we have explicitly written down 
the relation of the shear stress to the density perturbation 
in the synchronous gauge, after 5k in terms of 5c from 
Eq.(|21. 

The stress appears during the Cardassian take-over 
period, and increases temporarily the magnitude of the 
metric perturbation When the Cardassian density is 
beginning to dominate, the effect of stress is fading away. 
The integrated Sachs- Wolfe effect is stronger than in case 
I since now the effective matter source of the gravitational 
potentials is Pk^k instead of Pc5c- The dip in the poten- 
tial 4* may somewhat cancel the earlier contribution to 
the integrated Sachs- Wolfe effect, but the net effect is, as 
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FIG. 11: The total matter power spectra for different param- 
eter values in the case III. 



in case II, that the ISW again increases with increasing 
q. We show the evolution of gravitational potentials in 
FigElfor parameters (g, z^) = (1.5, 1). The CMB spectra 
is plotted in Fig^|for the same parameter choices as we 
did in case II. Choosing parameter values v ^ 1 seems to 
give similar results. 

We plot the total matter power spectrum in Figllll 
Evolution of the linear overdensities seems to be similar 
as in the Case I. Thus the shape of the matter power 
spectrum is not as sensitive to the parameters q and as 
the CMB spectrum at large scales. 



Conclusions 

We have attempted to describe modifications of the 
Friedmann equations in terms of a cosmological fluid with 



a priori unknown properties. We applied different as- 
sumptions on such a fluid to the representative case of 
the Cardassian model. In case I we assume that the fluid 
has no fluctuations and just provides a background for 
the evolution of the standard cosmological energy com- 
ponents. 

A more physical assumption is considered in case II 
where we interpret the Cardassian energy density as that 
of an interacting dark matter. It is shown that then 
the requirement of a physically tolerable sound speed re- 
stricts the model to the unmodified polytropic Cardas- 
sian one, which is equivalent to the generalized Chaplygin 
gas model. These are known to predict damping of the 
matter power spectrumj2Q| and CMB spectrum inconsis- 
tent with observations ^4j. Following ideas introduced 
in the Chaplygin gas literature (2"^.|28j|. we comment on 
the possibility to improve the results by modifying the 
fluid interpretation (of either Chaplygin gas scenario or 
the Cardassian expansion). 

In case III we interpret the additional energy associ- 
ated with the dark matter as arising from a modification 
to the general relativity. This is no longer compatible 
with perfectness of the fluid represented by T^. When 
Einstein equations are modified, we do not in general ex- 
pect equality of the Bardeen potentials in the absence of 
anisotropic stress. If we then write the Einstein equa- 
tions in the conventional form and take into account the 
modifications in the form of a fluid, it is generally found 
to be imperfect. This is also consistent with|29|, where 
modified Friedmann equations are investigated under the 
assumption that the Birkhoff theorem is satisfied. Then 
also vj/ and <& are found to differ. In fact our approach 
to modified gravity in the case III coincides with the one 
in 29] at the scales they consider. Note that in the formu- 
lation of our case III, there is no problem of interpreta- 
tion of the fluctuations in the present Universe. For the 
scales of interest, the dark matter perturbations evalu- 
ated in different gauges do coincide at late times. This 
is not true in the case II, as anticipated in^. Large 
scale structure in the limit fc — > in these models with 
assumptions equivalent to our case III was studied also 
in 30]. 

In both cases we find that the late integrated Sachs- 
Wolfe effect is typically very large in the Cardassian mod- 
els. For the fluid interpretation (case II) it is clear, also 
from the matter power spectrum, that the model is not 
compatible with observations unless one chooses param- 
eters very close to the ACDM model. In the description 
of modified gravity (case III) the matter power spectrum 
is not as drastically distorted when g, ^ 1, but it seems 
that the ISW effect could be used to rule out most of the 
parameter space also in this case. However, a detailed 
likelihood analysis is beyond the scope of the present pa- 
per. 

Hannestad and Mersini-Houghton^] have recently 
presented results of a more general investigation on the 
effects of new physics on the CMB spectrum. They also 
find that modifications of the late-time gravitational in- 
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teractions can lead to a similar boost of the low-^ part 
of the spectrum as seen here in the Cardassian model for 
our cases II and III. It remains to be seen how generic 
this effect is and what constraints it puts on the equation 
of state for the dark energy. 
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